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Abstract. Autonomous circadian clocks drive daily rhythms in physiology and behaviour. A network of 
coupled neurons, the suprachiasmatic nucleus (SCN), serves as a robust self-sustained circadian pacemaker. 



vN Synchronization of this timer to the environmental light-dark cycle is crucial for an organism's fitness. In 



o 



a recent theoretical and experimental study it was shown that coupling governs the entrainment range of 

"^^ circadian clocks. We apply the theory of coupled oscillators to analyse how diffusive and mean-field coupling 

d 

. ^ affects the entrainment range of interacting cells. Mean-field coupling leads to amplitude expansion of 

I I weak oscillators and, as a result, reduces the entrainment range. We also show that coupling determines 

T-H the rigidity of the synchronized SCN network, i.e. the relaxation rates upon perturbation. Our simulations 
> 

C and analytical calculations using generic oscillator models help to elucidate how coupling determines the 



in 

o 



X 



entrainment of the SCN. Our theoretical framework helps to interpret experimental data. 

PACS. XX.XX.XX No PACS code given 

PACS numbers: 05.45.Xt, 05. 45. -a, 87.18.Yt cleus located in the hypothalamus, the suprachiasmatic 

nucleus (SCN). SCN cells synchronize to each other and 

- , , . generate a robust ~ 24 h self-sustained oscillation that 

1 Introduction 

drives locomotor and hormonal daily rhythms even in the 
Circadian clocks are endogenous oscillators driving daily absence of external forcing. The environmental day-night 
rhythms in physiology and behaviour. In mammals, these cycles periodically modulate these oscillations [J. Thus 
rhythms are centrally controlled by a tiny neuronal nu- the circadian system can be regarded as a network of 



- We thank Werner Ebehng and Vadim Anishchenko for their ^^-^^^ ^^^^^^^ oscillators. Many details of the entrainment 
guidance through the world of coupled oscillators. 
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(synchronization to external environment) and intercellu- anisms are believed to be responsible for the SCN syn- 

lar synchronization processes are not well understood. In chronization: synaptic connections, gap junctions and se- 

particular, the coupling mechanisms between neurons are crated neuropeptides [10] . Quantitative details of the cou- 

debated and the coupling strengths are unknown. Here we pling mechanisms are still unknown |llj . In this work 

apply the established theory of coupled oscillators [HISlIll we study two basic types of coupling mechanisms: diffu- 

\5\ to study the synchronization and entrainment proper- sive and mean-field coupling. Diffusive coupling is gener- 

ties of the mammalian circadian clock. Our goal is a bet- ated by a difference in state between the given cell and 

ter understanding how the interplay between cell-to-cell its neighbourhood, hence its name, whereas in mean-field 

synchronization affects the entrainment properties of the coupling all cells equally contribute to the common mean- 



whole SCN network. field, which acts back on each of them (see Section 2.2 1. 

This paper is organized as follows: first we summarize 
the major relevant physiological findings on the entrain- 1.2 Entrainment of the circadian clock 

ment properties of the mammalian SCN and peripheral 

The circadian clock has evolved to synchronize an organ- 
tissues. In order to capture the most fundamental oscil- 

ism to periodically recurring environmental conditions, 
latory features of the circadian clock, we compare those 

such as light-dark or temperature cycles. Robust entrain- 
experimental findings with modeling results. We study for 

ment to external periodic signals is essential for a pre- 
this purpose rather generic coupled oscillators. We con- 
cise timing of behaviour and metabolism. External stimuli 
elude that observed differences in the entrainment range 

such as light or neuropeptide pulses can shift the phase of 
of SCN and peripheral tissues can be assigned to coupling- 

the SCN clock by a few hours [T2|[T3]. Periodically applied 
induced amplitude and rigidity changes. We discuss impli- 

external stimuli can entrain the SCN to a diverse range of 
cations for the interpretation of experimental results. 

periods typically within a range of 24 ± 2 h P^[T51[TB] . In 

addition to the SCN, almost every cell in the human body 
1.1 The SCN - a network of coupled of oscillators 

shows circadian oscillations, such oscillators are known as 

Oscillations in the SCN emerge at the individual cell level peripheral circadian oscillators |^[T7]. 
and are generated by intracellular genetic feedback loops In a recent series of experiments, temperature cycles 

[6]. Experiments with dispersed SCN neurons revealed within the physiological range were applied to organotypic 

that individual cells display oscillations with periods rang- SCN and lung slices and their entrainment ranges were 

ing from 20 to 28 hours 7 . Coupling between Individ- characterized [T51I19J . Even though the molecular mech- 

ual SCN neurons confers a precise and robust synchro- anism driving the oscillations at single cell level in pe- 

nized rhythm [8l[9]. Diverse intercellular coupling mech- ripheral and SCN tissues are quite similar [9lll7j. the lung 
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tissue was able to entrain to a much wider range |;19, . 
These different entrainment ranges reflect presumably dif- 
ferences in the intercellular coupling in the SCN and pe- 
ripheral tissues. Puzzled by these differences between the 
SCN and peripheral tissue clock properties, we analyze 
systematically how oscillator properties and coupling be- 
tween the oscillators affect the entrainment range. 

2 Results 

In the current paper, we make use of both direct numer- 
ical simulations of the system on hand and of numerical 
bifurcation analysis. Please see Appendix [C] for the details 
of the numerical methods. 

2.1 What oscillator properties govern the entrainment 
range? 

We follow the tradition of Winfree [5U], Kronauer [2T], 
and Glass & Mackey [2^ and study generic amplitude- 
phase oscillators. More detailed biochemical models p51 
I241I25J can be approximated by amplitude-phase models 
[2S] . Pure phase models [S] might be not sufhcient since it 
has been shown experimentally that the amplitude in the 
clock oscillations is variable [27ll28] and affects entrain- 
ment behaviour |29 l. l30 l [3T 1 [32] . 

A simple amplitude-phase oscillator can be generically 
written as 
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Fig. 1. Entrainment range for three different single oscillator 
models. In each frame, the borders of the entrainment range 
for rigid (7 — 1.0, black lines) and weak (7 — 10~^, red lines) 
oscillators with Ao = 1.0 and r = 24 are shown, (a) Oscillator 
with linear radial dynamics /; , (b) the Poincare oscillator, and 
(c) the Hopf oscillator. In all three frames (a), (b), and (c), 
the black lines occlude the red ones in the region with small 
b. Below red dots, red lines denote saddle-node bifurcation of 
limit cycles. Above red dots, red lines denote Neimark-Sacker 
bifurcation of limit cycles. Black lines always denote saddle- 
node bifurcation of limit cycles. 
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Here, the parameter r denotes the intrinsic period of the Figure [l] shows the entrainment ranges of amphtude- 

oscillator. The function / (r) determines the particular phase models described by Eq. S driven by a periodic 

type of oscillator. In this paper we consider three choices force of a period T: 

of/(r): 

z = (iuj^ -if (r)) z + 6e'^* (4) 

f ('') = fp (^) := ?■ — ^Oi ^ ' The parameter b denotes the Zeitgeber strength and T = 
f (r) — f] frl •= 1 — Anir 27r/i? denotes the Zeitgeberperiod. Depending on the Zeit- 
The nonhnearity A corresponds to a Hopf oscillator (which ^eber strength 6, the oscillator can synchronize to a range 
is sometimes called modified van der Pol oscillator, see ^^ Zeitgeber periods T, also known as range of entrain- 
[33]), fp corresponds to a Poincare oscillator [22], and fi °^^n^- Entrainment range becomes broader with increas- 
represents a linear dynamics in the radial variable. In all "^S ^- ^^us, we find the characteristic triangular shape of 
three cases, A^ determines the size of the hmit cycle. Pa- ^^'^ entrainment region, also termed 1:1 Arnold tongue [5]. 
rameter 7 determines the relaxation rate towards the limit ^or all three oscillator types, the entrainment range of the 
cycle with r = Ao. Equation Q can be written in a more ^^ak oscillator is broader than that of the rigid one, corn- 
compact complex form as P^^^ '^^'^ li"^^ ^S^i^^* black ones in Figure [l] 

We have identified the following mechanism that leads 

z ^ (iuj — ^ f iv]) z \ ) 

to the differences in the width of the entrainment range: 

for complex variable z = re"'' with amplitude r and phase Generically, the entrainment range of limit-cycle oscilla- 
ip. Here, the frequency of the oscillator is given hy uj = tors is bracketed by a pair of either saddle-node (SN) or 
2tt/t. In this subsection, the models refer to the SCN Neimark-Sacker (NS) bifurcation lines. In rigid oscillators, 
rhythm as a whole. Coupling of single cell oscillators will SN bifurcation lines (black lines in Figure [ij) continue up 
be discussed below. to relatively large values of & ^ 0.05, where the entrain- 
In the context of genetic circadian oscillators, we study ment range is about 10 hours wide. Contrarily, in weak 
mainly so-called weak oscillators, characterized by small oscillators, SN bifurcation is found only for small values 
positive 7, see a detailed quantification of clock oscillations of b (parts of red lines below red dots in Figure IT]) . With 
in [28]. As a consequence of small 7, the amplitude of weak b that small, SN lines of week and rigid oscillators coin- 
oscillators strongly depends on the applied forcing and/or cide, compare Figure IT] (b) and (c) with black lines De- 
coupling to neighbouring oscillators. As we will see, the eluding red ones below red dots. For larger values of b, 
changes in the oscillation amplitudes can cause shrinkage the entrainment range of weak oscillators is limited by NS 
or broadening of the entrainment range. lines (parts of red lines above red dots in Figure [T]), which 
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span a broader range in Zeitgeber period T. Thus, we at- 
tribute the broader entrainment range in weak oscillators 
to switching from SN bifurcation to NS bifurcation, which 
in our case is achieved for small 7 (rigid vs. weak oscilla- 
tors). Note that in general, the switching between SN to 
NS bifurcation can be controlled by different parameters. 

Saddle-node and Ncimark-Sacker bifurcation lines can — ■_'^ _■, ^ ■, 

cross in codimension-two Bogdanov-Takens points " those PGilOQ QITTGrGrlCG 1 T 




are exactly the red dots in Figure [T] An in-depth dis- 
cussion of such high-codimcnsion bifurcation points falls 
beyond the scope of the present paper and we refer to a 
standard textbook [5]. 



Fig. 2. Dependence of the amplitude of the entrained os- 
cillator on the detuning T — t. The solid red line refers to 
a rigid linear oscillator (7 = 1), the dashed black line repre- 
sents a weak linear oscillator (7 — 10"^), and the dotted black 

The amplitude of entrained oscillators within the en- ^'"^^ ■^e^^'^^ ^^ ^^ak Hopf oscillator (7 = 10'^). In each case, 

, . ; . 1 -1 u • TT rr>l A the forcing strength b is chosen in such a way that the width 

tramment range is exemplariiy shown m b igure ^ As we 00 j 



can see, weak oscillators demonstrate a pronounced am- 
plitude increase close to the middle of the entrainment 
range. Close to the limits of entrainment, the amplitude 
of weak oscillators drops below the value of the amplitude 
of the unperturbed oscillator [341. Amplitude changes in 
rigid oscillators do not show reasonable changes over the 
whole entrainment range. 



of the synchronization range is 4 h: for weak Hopf oscillator 
6 = 1.6 X 10"^, for rigid linear b = 2.2 x 10~^ and for weak 
linear & = 9 x 10"'^. The Poincare oscillator shows comparable 
behaviour. 



plitude Aq, and the relaxation rate 7) shape the entrain- 
ment range. 

In the following sections we take into account that the 
As intuitively expected, the ratio of Zeitgeber strength SCN oscillator is in fact a network of coupled single os- 
b to oscillator amplitude Aq determines the entrainment dilators. Consequently, we address the question how cou- 
range. This has been confirmed in |19| using extensive sim- pling influences the amplitude and the Floquet exponents 
ulations and analytical calculations. Figures [T] and |2] sug- of the synchronized SCN. We show that mean-field con- 
gest that the change in the entrainment range can be also pling reduces the entrainment range via amplitude expan- 
influenced by the entrained amplitude and the rigidity, sion. Moreover, for identical oscillators diffusive coupling 
Thus we conclude that for a given entrainment strength, affects rigidity but causes little effects on the amplitude 
all three oscillator properties (the intrinsic period r, am- and therefore no major changes in the entrainment range. 
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2.2 Mean-field coupling leads to amplitude expansion 

So far, we have considered the SCN as a single limit cycle 
oscillator. In the following we take into account that the 
SCN is actually a network of coupled cells. Under normal 
condition these cells are well synchronized. Specific condi- 
tions such as constant light [5S] or exotic short light-dark 
cycles might lead to desynchronization [36| . Here, we focus 
on synchronized cells driven by external rhythms. As dis- 
cussed above, multiple coupling mechanisms contribute to 
synchronization. We study here diffusive coupling, which 
represents, e.g., gap junctions, and mean- field coupling, 
which might model secreted neuropeptides such as VIP 

m- 

In case of two coupled cells these coupling mechanisms 
are described by the following equations for complex am- 
plitudes zi^2 = ri_2e"''i'^ of both oscillators 

ii = (iwi - 7/ (ri)) zi + ^,nf (zi + Z2) + ^J.d (z2 - zi) , 

Z2 = {1UJ2 - 7/ (^2)) Z2 + Mmf {Zl + Z2) + Hd {Zl - Z2) ■ 

(5) 

Here, ^1.2 = ^^^ are the oscillator frequencies with ri^2 
close to 24 h being their internal periods, ^nj accounts 
for the strength of mean-field coupling, and /id accounts 
for the strength of diffusive coupling. In this formulation, 
parameters /imf and /id can be varied independently. Pure 
diffusive coupling corresponds to /i,nf = and pure mean- 
field coupling to /id = 0. 

We have discussed in the previous section that the ra- 
tio of external forcing to oscillator amplitude determines 
the entrainment range. Thus for a given Zeii(/e6er strength 
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Fig. 3. (a) Color-coded standard deviation of phase dif- 
ference between two mean-field coupled Poincare oscillators 
(7 = 10~^,Ao = 1.0, Ti = 24). (b) Color-coded amplitude 
of two mean-field coupled Poincare oscillators. Amplitudes of 
both oscillator are identical. 

6, the amplitudes of the coupled oscillators are essential. 
Large amplitudes oscillator are difficult to entrain [2^[5T] . 

In identical oscillators that are completely synchro- 
nized (i.e. if Zl = Z2), the diffusive coupling term /id (zi — Z2) 
in Eq. ([5]) vanishes and, hence, the synchronized state is 
the same as in an uncoupled system. In other situations 
such as detuning of the frequencies, amplitude reduction 
can be achieved [3711381. 



Contrarily, pure mean-field coupling (i.e. /id = 0) can 
induce pronounced resonance behaviour. The coupling term 
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Mmf {zi + Z2) in Eq. (Is]) acts like a periodic driving force 2.3 Coupling strength determines rigidity of the 

and can enhance in this way oscillation amplitudes of zi coupled system 
and Z2 drastically. This effect is particularly strong if the 



We have seen in section [2?T] that beside period and ampli- 
oscillators are weak (small 7) and especially for Poincare 

tude, relaxation rates (the rigidity of oscillators) affect the 
oscillator and oscillator with linearized dynamics, compare 

entrainment range. Below we will demonstrate that in a 
section 12.11 

network of coupled oscillators the slowest relaxation rates 

(the Floquet exponents of the synchronized system as a 

We derive in the Appendix A an analytical expression 

whole) are governed by the coupling strength. 

for the amplitude of the synchronized system. It turns out 

Let us consider formally an ensemble of N identical 

that mean-field coupling generally expands amplitudes, ,1 .1, mi ■ 1 1 • 1 

uncoupled oscillators. This large system obviously pos- 

particularly in the middle of the Arnold tongue. Contrar- . , . . , 

sesses N zero Floquet exponents which are just the trivial 

ily, diffusive coupling reduces the amplitude and might 

Floquet exponents of the member oscillators. Now, by in- 
even suspend the oscillations completely. Both effects are 

troducing a small coupling, we expect that the N zero 

pronounced for weak oscillators, i.e. for small values of 7. r 1 i t i i 

b loquet exponents of the huge system slightly move away 

from zero, while remaining in a small vicinity of zero (ex- 
Such an amplitude expansion under mean-field cou- ^^^^ ^^ ^^g ^ero exponent which remains at zero due to 
phng is illustrated in Figure [3](b). The 1:1 Arnold tongue ^he phase shift invariance). For smaU couphng, the shift 
is clearly visible in the plot of the phase difference in Fig- ^f Floquet exponents from zero will be proportional to 
ure [3] (a). Figure [s] (b) shows an amplitude expansion by ^^e coupling strength. Thus, the slowest time scale in the 
two orders of magnitude, particularly within the Arnold system will be dictated by the interaction between oscil- 



tongue. In section 2.4 we will exemplify how such an am- lators 

plitude expansion leads to a drastic shrinkage of the en- j^ Appendix B we present the matrix M whose eigen- 

tramment range. values approximate the Floquet exponents of the synchro- 
nized state of weakly coupled oscillators. This matrix has 

In the following section we will demonstrate that cou- the following properties: It applies to a general setup of an 

pling strength is also intimately related to the rigidity of arbitrary number of limit-cycle oscillators of any dimen- 

the coupled system as a whole. Consequently, coupling sion and is universally applicable to different models of 

might control the entrainment range in two ways: Via circadian clocks, including the simple models considered 

amplitude (the results of this section) and via rigidity as in this paper. Matrix M can be deduced just from the 

discussed below. properties of the unperturbed limit cycle oscillator. This 
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in turn implies that using matrix Af , it is possible to cal- size of the limit cycle and thus the increasing limit cycle 

culate the rigidity of the coupled system as a whole only becomes more stable. This is of course not the case for 

from information on the single uncoupled oscillator. This diffusive coupling, since diffusive coupling does not lead 

matrix depends only on how the coupled oscillators cross- to an amplitude increase (see Figure l4] (a)), 
talk to each other. In this sense, this matrix is applicable 
for both mean-field and diffusive coupling. 

As a main result, we find that the relaxation rates 2.4 Coupling governs entrainment range 

(the Floquet exponents) of the synchronized system are 
proportional to the coupling strength for almost any os- 
cillator model and coupling type. We have shown above that coupling strength determines 

„. m-ii . . .1 T 1 rm , amplitudes and Floquet exponents of a coupled oscillator 
h igure |4| illustrates the Imear decay of t loquet expo- 

. , .. , x^ 1-rr • T • system. It has been argued in Introduction that entrain- 
nents with coupling strength, i'or diffusive coupling as m 

„. m / \ 1 , 1 • 1 -,1 ,, ment to external rhythms is crucial for an organism's fit- 
I'lgure WJ (a), homogeneous perturbations decay with the 

,,, , . ,11. . TT rr\/ \\ ■ ,1 ness. It was found experimentally that coupled SCN neu- 
ratc 7 (the horizontal line m l'igureWJ(a)j, since the cou- 

,. • , 1 , A rons are harder to entrain than peripherial tissues [T51I19I . 
phng terms 2:1^2 — ^2,1 vanish due to the symmetry. As 

.,.,., ; 1 1 , , 1 , • 1 We illustrate in this section that rigidity and amplitude 
intuitively expected, heterogeneous perturbations decay 

.,, , ,. , , ;, ,. , ,1 ^ expansion affect the widths of the Arnold tongues, 
with a rate proportional to the coupling strength. Conse- 

quently, monitoring of transients in experiments can po- Figure [5] shows the entrainment ranges of coupled rigid 
tentially provide information on coupling strength. Fig- (a, b) and weak (c,d) oscillators. As a model, we took cou- 
ure I4] (b) shows the Floquet exponents of the mean-field pled Eqs. (l5| with right-hand sides for zi and Z2 perturbed 
coupled system. Beside the trivial exponent at zero, there by the same periodic forcing in the form of 6e'T^* with 
are two small negative exponents that grow into negative the amplitude b and period T. It turns out that coupled 
values as the coupling strength /imf increases. In contrast rigid cells in Figure p^ (b) show no amplitude expansion 
to Figure [4] (a), the single oscillator relaxation rate 7 does and, consequently, the Arnold tongue resembles the en- 
not persist at a constant value as coupling strength in- trainment zones in Figure IT] In contrast, amplitude ex- 
creases. We explain this growth of the absolute value of pansion of weak oscillators leads to a drastic shrinkage of 
the exponents by the fact that increasing mean- field cou- the Arnold tongue, see Figure Isjc) and (d). Similar effects 
pling leads to an increase of the amplitude of the syn- of mean-field coupling are observed in the other oscillator 
chronized state (see for instance Figure [31(b)). Due to the models (the Hopf oscillator and the oscillator with linear 
nonlinearity of the system, the exponent depends on the radial dynamics). 
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3 Discussion 

Circadian oscillators have strikingly similar molecular mech- 
anisms both in SCN neurons and in peripheral cells [51155] . 
At single cell level, the relative amplitudes vary from 0.1 
to 2 and estimated damping rates are in the range from 

. Single cell periods have a stan- 
dard deviation of 1-2 h [71[S]. Despite such variable and 
noisy single cell oscillations, the SCN as a whole is a re- 
ally precise pacemaker with a day-to-day variation of a 
few minutes [S]. This precision is believed to be achieved 
by intercellular coupling of SCN neurons jTT]. 

Recent experiments revealed that coupling governs the 
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0.00 0.02 0.04 0.06 0.08 0.1 entrainment range of circadian clocks [T51IT5] . Lung tis- 

sue without strong intercellular coupling was entrained 
Fig. 4. Dependence of the Floquet exponents A on coupling by 1.5°C temperature cycles at external periods of 20 h 
strength, (a) Two diffusively coupled (strenth ^d) Poincare os- and 28 h. Contrarily, SCN tissue could not be entrained 
dilators with 7 = 10~^. Dashed horizontal line denotes the even by a larger temperature amplitude of 4°C and an 
local Floquet exponent which is not affected by the coupling external period of 22 h 19J. Only if the coupling was re- 
strength. In contrast, the other sloped dashed line denote the duced via chemicals (MDL and TTX), entrainment was 
interaction Floquet exponent that decreases with increasing achieved. These experimental results motivated our sys- 
couphng strength, (b) Two mean-field coupled (strength ^mf) tematic analysis of coupling and entrainment. 



Hopf oscillators with 7 = 10^^. In both (a) and (b), red 
straight lines denote the trivial Floquet exponent at zero. Due 
to the amplitude expansions, both non-trivial Floquet multi- 
pliers (dashed lines) decrease as coupling strength increases. 
Also, there always exists a fourth exponent at much smaller 
values of A, which is not included in (a) and (b). 



We have shown in this paper that coupling can affect 
the oscillator properties of the network in two ways: (i) via 
amplitude expansion due to mean-field coupling and (ii) 
via increased rigidity of the SCN network. Both effects 
have been suggested to explain the differences between 
lung and SCN tissues [T5] . 

Our results on the Floquet multipliers of the synchro- 
nized state in ensembles of oscillators have been obtained 
in a quite general setup for an arbitrary number of oscil- 
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lators and connection topology. We also strongly believe are, for instance, indications that coupled SCN neurons 

that our numerical results can be straight-forwardly gen- exhibit larger amplitudes than dissociated neurons 0[1S] . 

eralized to a larger number of interacting oscillators. This could reflect amplitude expansion via mean-field cou- 

On the other hand, the structure of Arnold tongue be- phng as studied in section [2. 2[ The enlarged entrainment 

comes more involved for oscillator ensembles with different range due to MDL treatment in 19 might be related to 

internal frequencies. Figure [5] suggest that already two os- reduced rigidity, since the relative amplitudes are unaf- 

cillators result in a W-shaped entrainment range. We also fected. 

speculate that the visible rippling in Figure [5] might be In summary, our results suggest that the established 

attributed to a secondary bifurcation structure due to the theory of coupled oscillators can provide valuable insights 

mismatch of the internal periods. in the field of circadian rhythms. 

For a quantitative comparison of our simulations with 
the in vivo SCN network, direct measurements of the cou- 

A Entrained amplitude for mean-field and 

pling strength and oscillator rigidity are desirable. How- 

., , , rr., 1 diffusive coupling 

ever, such data are currently not available. There are chem- 
ical treatments with TTX [ST and MDL gO], studies with 

The aim of this appendix is to demonstrate how mean- 
dissociated neurons [5112] and mechanically cut slices [27 . 

field and diffusional coupling together with the amount of 

Furthermore, knockout studies have been performed 

detuning affects the amplitude of the synchronized solu- 

extensively (VIP, VIPR, gap junctions). Unfortunately, all 

tion. Our calculations are based on the method described 

these interventions lead to considerable destruction of the 

in [38]. 

SCN network and to poorly controlled side effects such as 

Rewriting Eq. (ISJ) in polar form with zi.2 = ri^2^^f^-'^ 

changes of ion and neurotransmitter concentrations. 

and looking for a stationary solution with ri = r2 = r, we 

Thus coupling mechanisms in the intact SCN have to 

arrive at 

be explored indirectly. Using powerful imaging techniques, 

r = (Mm + A^p cos 9 -"ff (r)) r, 
amplitudes and phases of single cells can be monitored [HI (6) 

— — — ini ^ • I- • ^ /A — 2/ir) SlU P, 

119] . lemperature cycles allow the analysis of entrainment 

properties. Moreover, temperature pulse response can be where we have defined 9 ^ ipi - ip2, A = loi - LO2, ^,p = 

used to derive phase response curves [151141] . Such indirect /^mf + Md, and /i„j — ^i-^^i — /id. 

measurements can be exploited to infer properties of single Looking for phase- locking regimes is equivalent io 9 = 

cell oscillators and their coupling. This reverse approach 0, which immediately results in sin0 = -^. We can use 

can profit from the theory presented in this paper. There this to express cos 9 in the first equation of Eq. ([6|) and 
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Fig. 5. Phase difference and amplitude behaviour of two mean-field coupled Poincare oscillators under external forcing. Color- 
coded standard deviation of the phase difference (left panels) and entrained amplitude (right panels) . Both coupled oscillators 
are subjected to the same external Zeitgeber period T and Zeitgeber strength h. Further parameters are: A012 = l,Mmf ~ 
0.005, Ti — 23.5, T2 — 24.5, (a) and (b) rigid oscillators (71,2 ~ !)• Note the W-shaped Arnold tongue due to the fact that 
ri 7^ T2 and hence the individual entrainment ranges of both oscillators are shifted with respect to each other, (c) and (d) weak 
oscillators (71,2 ~ 10"'^). 



impose f = 0, which results in either r = or 



/(r)=Um±VM^-T /7 



one of the solution (with plus sign) is a stable one, whereas 
the second one is unstable. 

(7) Some special cases are worth mentioning here. Con- 

sider first two identical oscillators with A = 0, then the 
above expression for the amplitude r simplifies to 



/ ('') = (Mm ± Aip) /l 



Here, / (r) is one of the three possible radial dynamics 

from Eq. (l2|. Thus, Eq. (It]) constitutes the conditional 

equation for the amplitude r of the synchronized state. 

We see that both coefficients for mean-field and diffusive which results in either / (r) = 2/i,jjf/7 or / (r) = 2^^/j. 

coupling contribute to the amplitude change. Note that This explains the amplitude expansion for the case of pure 
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mean-field coupling with /id ~ 0. Moreover, choosing yim,p calculation in small parameter e shows that the N y. (.- 

such that (/im + jJ-p) ll — ~1 "kills" the linearly unsta- dimensional limit cycle oscillator exhibits A^ — 1 Floquet 

ble part of / (r) and thus results in an absence of oscilla- eigenvalues being e-small. One Floquet eigenvalue is zero, 

tions |37j . as implied by the phase shift invariance. The N critical 

The condition for a phase-locking solution sin = ^ (those close to zero) Floquet exponents of Eq. Q are ap- 

describes the width of the Arnold tongue, giving — 2/ip < proximated to the e order by the eigenvalues of the foUow- 

A < 2fip. At the border of the entrainment range, both ing matrix M: 
stable and unstable solutions disappear in a saddle-node 



bifurcation. For given fj,p and firm the amplitude of the 
stable solution (with plus sign) has a maximum at Z\ = 0, 
which follows from the inspection of Eq. (ItI) . ]\/[ 

B Floquet exponents of fully synchronized 
state of identical oscillators 



i 



Kjvi Kjv2 ... — 2_^ '^Ni 



(10) 



Let us consider N identical oscillators, each of them de- 
scribed by a ^-dimensional state vector Ui e M^, i — 
1, 2, . . . , iV. Suppose that in the absence of coupling each 
of the oscillators is described by the equation 



The coefficients in the matrix depend on the integrals 



dtij{t)K,j{t-e)q{t), 





where ij: (t) is the eigenfunction to the zero eigenvalue 
u = g [u) . [°) of the adjoint linearization of Eq. ([8| around the unper- 



The unperturbed limit cycle is u{t) ~ q (t) with period r 
such that q{t + t) = q (t) for all t. Now we shall study a 
coupled ensemble of N identical oscillators 



turbed limit cycle q {t) and Kij — |^ are the Jacobian 
matrix elements of the coupling functions ki evaluated 
along q{t). Note that the matrix M does not contain 
terms k^. Those reflect the influence of the Ui on itself 
Ui= g (ui) + ki (ui, U2, . . . , UN', e) , i ~ 1,2,. . . ,N, (9) through the couphng function ki. Thus, for both mean- 
field coupling and diffusive coupling the matrix M will be 

where ki (ui, U2, . . . , un', e) are coupling functions depend- 

the same. 

ing on the coupling strength e. For small e, the coupling 

functions ki are also assumed to be small, i.e. fc^ (. . . ; e) = If all coupling functions ki {ui,U2, . . . , un', e) depend 

O (e). A straight-forward, but rather tedious perturbation on aU other Uj,j 7^ i in the same manner, that is, all g^ 
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are equal, the matrix M simplifies to 



/ 



M = K„ 



a 13 



-a 






13 /3 



with 



/dk 
dtxl:{t)-^{t-e)q{t), 



and a = {N -l)/N and j3 = l/N. For ki{ui, . . . ,UN;e) = 
eN-^Y.f MqUj we have |^ = eN-Hio for all i ^ j, 
where Mo is a i? x ^ matrix with constant coefficients that 
describe how each of £ components of u are coupled to 
each other. Hence, Kq is just a scalar proportional to e. 

The eigenvalues of M are given by zero and iV — 1 
eigenvalues at —Kq- The eigenvector to the zero eigen- 
value corresponds to a simultaneous shift of all N phases 
by the same amount. The eigenvectors to the non-zero 
eigenvalues are given by the periodic shifts of the vector 
{N - 1, -1,-1, . . . , -1)"^. In the limit of large N those 
eigenvectors correspond to the "evaporation" eigenvalue: 
they describe the stability of the synchronized state against 
perturbations of a single oscillator, see |42j . 

C Numerical Methods 

Direct numerical simulations were performed with the help 
of the explicit 4th-order Runge-Kutta method with a time 
step of At ~ 0.01 d.u. (compare that to the typical period 
of oscillation t — 24 d.u.). Averaging for each parameter 
point in Figures |3] and [5] was done over the last 40 periods 
of total 140 periods of simulated oscillations. Bifurcation 



lines in Figures [T] [2J and Floquet exponents in Figure 
[4] were obtained with the help of the standard numeri- 
cal continuation and bifurcation analysis software AUTO 

2000 gg. 
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